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As a probe of the Yang-Mills vacuum, we study numerically the eigenmode spectrum of the covariant lattice 
Laplacian operator. We find that the eigenmodes at the low and high ends of the spectrum are localized in 
finite regions whose volume is insensitive to the lattice volume. We also find that the vacuum is seen very 
differently by localized modes of the covariant Laplacian in different representations of the gauge group. In 
the fundamental representation, the data suggests that the localization volume is finite in physical units set by 
the string tension, and localization disappears when center vortices are removed. In the adjoint and j = 3/2 
representations the low and high-lying modes are far more localized, and the localization volume appears to 
scale to zero, in physical units, in the continuum limit. The adjoint Laplacian is insensitive to vortex removal, 
but we find that exponential localization is absent for adjoint eigenmodes in the Higgs phase of a gauge-Higgs 
theory. Localization is also absent in the spectrum of the Coulomb gauge Faddeev-Popov operator, as required 
in Coulomb gauge confinement scenarios. 



I. INTRODUCTION 

The zero and low-lying eigenmodes of the lattice Dirac op- 
erator are of considerable interest in lattice QCD, because of 
the connection of low-lying modes (via the Banks-Casher re- 
lation [1]) to the chiral condensate, and the connection of zero 
modes (via the Atiyah-Singer index theorem [2]) to topolog- 
ical charge. It has been shown recently by Golterman and 
Shamir [3], for the two-flavor Wilson-Dirac operator outside 
the Aoki phase, that the low-lying eigenmodes are localized 
in finite regions which do not grow with lattice volume. The 
phenomenon is clearly related to Anderson localization [4] of 
particle energy eigenstates in a disordered medium. On the 
other hand, Aubin et al. [5] find that the zero modes of the 
Asqtad [6] lattice Dirac operator extend throughout the lat- 
tice, and yet concentrate in regions which scale with coupling 
like manifolds of effective dimension between two and three. 1 
It is conceivable that this finding relates to the suggestion [8] 
that thin brane-like objects must exist in the QCD vacuum, a 
suggestion which is supported by studies of the action density 
at P- vortex locations [9]. 

It would be interesting to know whether the localization (or 
concentration) of low-lying eigenmodes is peculiar to lattice 
Dirac operators, or whether the localization feature is generic 
in confining lattice gauge theories, and could perhaps serve 
as another type of order parameter for confinement. To ad- 
dress this question, we study in this article the localization 
properties of the spectrum of the covariant lattice Laplacian in 
the fundamental, adjoint, and j = 3/2 representations of the 
SU(2) gauge group. Unlike the Dirac operator, the covariant 



In this connection, we should also take note of the rather different results 
obtained by the Kentucky group [7], who compute topological charge den- 
sity via the overlap operator. Their finding is that most of the lattice volume 
is taken up by two space-filling folded sheets of opposite topological charge 
density, and dimensionality D = 3. 



Laplacian is straightforward to construct on the lattice, having 
neither chiral symmetries which ought to be respected, nor 
doublers which need to be removed. 

Our main finding will be that there are indeed localized 
eigenmodes of the covariant Laplacian at the extreme lower 
and upper ends of the spectrum, and we will argue below that 
localization implies that neither the physical mass of scalar 
particles, nor the mass of scalar particle bound states, can 
be adjusted to zero via the bare mass term, at least in the 
quenched approximation. However, a very unexpected result 
of our investigation concerns the qualitative dependence of lo- 
calization on the gauge group representation of the covariant 
Laplacian. In the fundamental j = j representation of color 
SU(2), localized states appear to concentrate in 4-volumes of 
finite extension in physical units, and this localization disap- 
pears when center vortices are removed from the gauge field 
configurations. In sharp contrast, in the adjoint and j = 3/2 
representation, the low-lying eigenmodes are localized in far 
smaller regions, whose volume seems to go to zero in physical 
units in the continuum limit. 

We also find that localization is absent in the low-lying 
eigenmodes of the Coulomb gauge Faddeev-Popov opera- 
tor, despite that operator's close resemblance to the covariant 
Laplacian. This absence of localization is implied by the exis- 
tence of a confining color Coulomb potential, and is therefore 
required in Coulomb gauge confinement scenarios. 

In section 2, below, we present our results for the Laplacian 
eigenmodes in the fundamental representation. Section 3 is 
concerned with the eigenmodes of the adjoint representation 
Laplacian, in both pure Yang-Mills and gauge-Higgs theory. 
Results for the j = 3/2 representation are presented in section 
4. Section 5 concerns the relation between localized modes 
and the absence of massless particle states; this section also 
contains our results for the Coulomb gauge Faddeev-Popov 
operator. Our conclusions are in section 6. 
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II. LOCALIZATION IN THE FUNDAMENTAL 
REPRESENTATION 

The covariant lattice Laplacian in the fundamental repre- 
sentation of the gauge group is given by 



-VV 



(2.1) 



(color indices a,b for SU(2) are 1-2) and we are interested 
in the low-lying eigenmodes §%{x) satisfying the eigenvalue 
equation 

-A$ti(y)=X n tt{x) (2.2) 

Both the eigenmodes and eigenvalues are, of course, depen- 
dent on the configuration U M (x). The observables relevant for 
localization are the Inverse Participation Ratios (IPR), defined 

as 



IPR n = V(£pl(x)) 



where p„(x) is the normalized eigenmode density 
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p«w=lcw 



(2.3) 



(2.4) 



with Pn{x) = 1, and V = Y, x 1 is the lattice volume in lattice 
units. If the n-th eigenmode is localized, with p(x) ~ l/b in 
some volume b and p(x) «0 elsewhere, then the IPR will 
grow linearly with lattice volume, on a hypercubic lattice of 
length L, as 



L 4 

IPRn = -r 
b 



(2.5) 



In contrast, if the eigenmode is extended over the full lattice 
with p(x) ~ 1 /L 4 , then the IPR will not grow with volume. 

Of course, intermediate situations in which the eigenmode 
is partially localized are also possible. It could be, for ex- 
ample, that while the eigenmode is largely concentrated in 
some region of volume b, still there is a substantial fraction 
of the norm which is accounted for by the bulk of the lat- 
tice volume. For this reason it is useful to study the falloff 
in p(x) away from its region of maximum concentration. We 
introduce the "Remaining Norm" (RN) observable, defined 
for a given eigenmode as follows: Place the values of p„(x) 
into a 1-dimensional array r(i), i = 1,2, ..,V, sorted so that 
r{i) > r(i + 1 ) . Then the remaining norm is given by 



R(K) 



(=1 



(2.6) 



The RN is the amount of the total norm (= 1) remaining after 
counting contributions from the subset of K < V sites having 
the largest values of p(x). If the curve R(K) becomes more 
extended as the lattice volume increases, without approaching 
some limiting envelope, then the eigenmode is not really lo- 
calized in a fixed volume around a single site. It should be 
noted that a closely related quantity, 



GW = LK0 



(2.7) 
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FIG. 1: Log-log plot of IPR — A vs. lattice length L, for the lowest 
eigenmode of the covariant Laplacian in the fundamental representa- 
tion, at P = 2. 1 - 2.4. 
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FIG. 2: Log-log plot of IPR— A vs. lattice volume in physical units 
(fm 4 ). The dotted line is a best fit to linear dependence on volume. 



where x = K/V is the volume fraction, was introduced previ- 
ously by Horvath in the second article of ref. [7], where it is 
called the "normalized cumulative function of intensity." 

We have calculated the average IPR of the lowest-lying 
eigenmode at {$ = 2. 1 , 2.2, 2 .3 , 2.4, and fit the data at each cou- 
pling to 



IPR(fi,L)=A(fi) + 



*(P) 



(2.8) 



i=\ 



where L is the lattice length in lattice units. By the "sub- 
tracted" IPR data we refer to our data for IPR($,L) minus 
the constant A(p) obtained from the best fit to eq. (2.8). Fig- 
ure 1 is a log-log plot of the subtracted IPR vs. L, together 
with the best fit. The proportionality of IPR— A(p) to lat- 
tice volume is quite clear from the match of data points to the 
fit. In Fig. 2 we plot the subtracted IPR data vs. the lattice 
volume V in physical units at a variety of couplings. Here 
V[fm 4 ] = Va 4 (p) = (MP)) 4 is the lattice volume in physi- 
cal units, where a — a(p) denotes the lattice spacing (in fm) at 
coupling p . In these units we find that the IPR values at differ- 
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ent p fall roughly on the same straight line, which implies that 
the localization 4-volume (ba 4 ) is constant in physical units. 2 
The next question is whether the localization we have found 
is due specifically to the confining disorder of the lattice con- 
figuration, or whether it is due simply to local stochastic 
fluctuations in the link variables. To address this question, 
we have computed the low-lying eigenmode spectrum in the 
center-projected, and vortex-removed configurations. Recall 
that the vortex-removed configurations [10] are obtained by 
multiplying the lattice configuration in maximal center gauge 
by the center-projected configuration. This is known to re- 
move the confinement properties of the lattice configuration 
(the asymptotic string tension is zero), but is otherwise a rather 
minimal disturbance of the original lattice, since plaquette val- 
ues away from P-vortex locations are left unchanged. Locally, 
the vortex-removed configuration is actually somewhat more 
disordered than the original configuration, since the plaque- 
tte values at P-vortex locations are multiplied by — 1 . Figure 
3 shows that the lowest eigenmode computed in the center- 
projected (or "vortex-only") configuration is localized much 
as in the original lattice. The vortex-removed data, however, 
is drastically different. In Fig. 4, we see that that the IPRs 
of the vortex-removed configurations remain essentially con- 
stant as L increases, and therefore localization is completely 
absent. The obvious conclusion is that localization is associ- 
ated with large-scale, confining disorder provided by center 
vortices, rather than local, non-confining fluctuations of the 
link variables. 
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FIG. 3: Subtracted IPR of the lowest Laplacian eigenmode vs. lat- 
tice volume V in physical units, for the "vortex-only" (i.e. center- 
projected) lattice configurations. 
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FIG. 4: Same as Fig. 3 for the configurations with vortices removed. 
In this case the data is consistent with extended, rather than localized, 
eigenmodes. 



ing Norm observable, defined in eq. (2.6). Figures 5-7 show 
the Remaining Norms at (3 = 2.1 and at various lattice vol- 
umes, for unmodified, center-projected, and vortex-removed 
configurations. For the unmodified configurations (Fig. 5), the 
curve does seem to become slightly broader as the lattice vol- 
ume increases, although the indications are that the RN data 
converges to a limiting curve at the largest volumes. The con- 
vergence is much clearer for the center-projected configura- 
tions, where almost all the RN data falls on essentially the 
same curve (Fig. 6). For the vortex-removed data, there is no 
convergence at all; the curve simply broadens as the volume 
increases (Fig. 7). 3 
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The same conclusion is reached by examining the Remain- 



FIG. 5: Remaining Norm vs. volume at (3 = 2.1 for unmodified lat- 
tices. 



2 In this figure we have included some IPR data at p > 2.4 obtained on only 
a single volume: L = 24 for P = 2.45,2.5, and L = 28 at p = 2.55,2.6. 
In these cases, which do not allow a fit to eq. (2.8), we have subtracted 
A(P) = 1.38 obtained at p = 2.4. Although these higher-p data points are 
consistent with the multi-volume data obtained for P < 2.4, they come with 
the caveat that the subtraction constant A(P) was not obtained from a multi- 
volume fit. 



Finally, in Fig. 8 we display the extension b l l 4 a of the local- 
ized region in physical units, for both the unmodified and the 
center-projected configurations, at couplings p = 2.1 —2.4. 



3 If we plot the vortex-removed data vs. volume fraction rather than volume, 
all the curves at different volumes coincide. 
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FIG. 6: Remaining Norm vs. volume at (3 = 2. 1 for center-projected 
lattices. 
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FIG. 7: Remaining Norm vs. volume at P = 2.1 for vortex-removed 
lattices. 



The localization volume b, in lattice units, is extracted from 
a best fit of the IPR data to eq. (2.8). It seems clear that the 
average extension of the localized region in physical units is 
essentially p-independent, with this length roughly 2.3 fm for 
unmodified configurations, and 1 .5 fm for vortex-only config- 
urations. In Fig. 9 we display, for individual lattice volumes 
at some higher P values, a quantity which would equal (ac- 
cording to (2.8)) the inverse localization volume (ba 4 )^ 1 in 
physical units. Again we observe from the data that this quan- 
tity is insensitive to the value of p (with the same caveats as 
in footnote 2 above). 

So far we have considered only the lowest eigenmode. For 
the higher modes, we find at any fixed volume that the IPR„ 
decreases (b increases) as the eigenvalue X n increases , until 
there is no localization at all. For the non-localized states, we 
see that the IPR are small (O(l)) even on large lattices, and 
do not increase with lattice size (Fig. 10). A similar result is 
obtained in the center-projected lattices (Fig. 11). This situa- 
tion is familiar in the Anderson model [4], where eigenmodes 
at the upper and lower ends of the spectrum are localized, and 
eigenmodes in the bulk of the spectrum are extended. The 
boundaries between localized and extended states are known 




FIG. 8: Extension of the localization volume in physical units. 
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FIG. 9: Inverse localization volume in physical units (see text). 



as "mobility edges" [11]. In our case we observe numerically 
a remarkable symmetry: If the eigenmodes (j)„ are ordered 
from lowest (n = 1) to highest (n = n max = 2L 4 ) eigenvalue, 
then we find for each lattice configuration U fl (x) the corre- 
spondence 



IPR„ - IPRn max -n+l 



(2.9) 



We emphasize that this identity holds for IPR's calculated in 
each configuration, and not just for the expectation values. 
The identity follows from a symmetry. Let (j)^(x) be an eigen- 
mode of the covariant Laplacian with eigenvalue 



X„=2D- <E n 



(2.10) 



where — £„ is the eigenvalue of the nearest-neighbor piece of 
the Laplacian. Then it is not hard to see that 



f = (-l) E "^„ 



(2.11) 



is also an eigenstate of the covariant Laplacian, with eigen- 
value 



X' = 2D + <E„ 



(2.12) 



Thus there is a pairwise correspondence of eigenstates with 
eigenvalues below and above X = 2D, having exactly the same 
IPR values. This accounts for the identity (2.9). 
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FIG. 10: IPRs of the first hundred eigenmodes of the covariant Lapla- 
cian operator at (3 = 2.1. 
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FIG. 11: Same as Fig. 10 for the center-projected configurations. 



III. LOCALIZATION IN THE ADJOINT 
REPRESENTATION 

One might expect that merely changing the gauge group 
representation of the covariant Laplacian operator, from fun- 
damental to adjoint, would not greatly affect the localization 
properties of the eigenmodes; certainly one does not expect to 
see any qualitative difference. Surprisingly, this expectation 
turns out to be wrong. 

The covariant Laplacian operator in the adjoint representa- 
tion of SU(2) has the same form in as eq. (2. 1), except that the 
color indices run from 1 to 3, and the link matrices 



Uf= l 1 Tr[a a U (l (x)o b U^x)} 



(3.1) 



are SO(3) group-valued (where the link variables inside the 
trace are in fundamental representation). The results for the 
IPR of the lowest eigenmode are shown in Figs. 12 and 15. 
The first thing to notice is that in comparison to the IPRs of the 
fundamental representation in Fig. 1, the corresponding IPRs 
of the adjoint representation are far higher. The difference, at 
fixed volume and fixed coupling p, is a factor of at least an 
order of magnitude, and this factor grows with p. Moreover, 



in contrast to the fundamental representation, the profile of 
the Remaining Norm actually shrinks with increasing volume 
(cf. Figs. 13 and 14 for p = 2.3), indicating that the lowest 
eigenmode becomes, if anything, more localized as the lattice 
volume increases. 4 
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FIG. 12: Log-log plot of the subtracted IPR vs. lattice length L, for 
the lowest eigenmode of the covariant Laplacian in the adjoint repre- 
sentation, at (3 = 2. 1 - 2.4. 
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FIG. 13: Remaining Norm for the lowest eigenmode of the covariant 
Laplacian, in the adjoint representation, for various lattice volumes 
at (3 = 2.3. Only the RN data for small numbers of sites are shown. 

However, if we extract Z?(p) from a fit to eq. (2.8), we find 
that ba 4 falls rapidly with P, indicating that the localization 
volume shrinks in physical units. Instead, if we compute ba 2 
at each P, it is found that this quantity is essentially constant. 
In Fig. 15 we have plotted IPR x a 2 versus lattice volume in 
physical units, and it is clear that the data points at different 
p-values fall on the same line. This would imply, from eq. 
(2.8), that ba 2 rather than ba 4 is p-independent, which can be 



4 The high degree of localization present in low-lying adjoint representation 
eigenmodes was pointed out several years ago by Bertie et al. [12], in con- 
nection with a Laplacian version of direct maximal center gauge. 
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Covariant adjoint Laplacian, p=2.3 
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FIG. 14: Remaining Norm for the lowest eigenmode of the covariant 
Laplacian, in the adjoint representation, for various lattice volumes 
at P = 2.3. All RN data are shown. 
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FIG. 15: Log-log plot of adjoint representation IPR x a 2 vs. lattice 
volume in physical units (fm 4 ). 



seen directly in Fig. 16. We also display, in Fig. 17, the data 
for IPR /(a 2 L 4 ) at some individual volumes at higher p. These 
points should approximate 1 / (ba 2 ), and again the data is only 
very weakly dependent on p. 

The fact that ba 2 is constant, rather than ba 4 , is open to 
interpretation. It could be, following the suggestion in ref. 
[5], that the lowest-lying eigenmode has its support on a 
surface-like region of area A p h ys , whose expectation value is 
p-independent, and lattice scale thickness. If the thickness 
is, for example, only a single lattice spacing, then in physical 
units 



From fits IPR(P)=A(P)+V/b(P), cov. adjoint Laplacian 



ba 4 = A phys a 2 



(3.2) 



Thus, if ba 2 is P-independent as p — > °°, that means that the 
surface area A p h ys has a finite extension in physical units. This 
in turn implies that the lowest eigenmode of the covariant 
Laplacian, in adjoint representation, has its support on a 2- 
brane of finite extension. Assuming minimal (one lattice spac- 
ing) thickness, one can estimate from O p i Jys — (440 MeV) 2 and 
Fig. 16 that the area of this brane is A p ^ ys s=a 26 fm 2 . This pos- 
sibility is attractive, because of the possible connection with 
ultraviolet renormalons [8] and/or P- vortex surfaces [9]. 

However, a more straightforward interpretation of Fig. 16 
is that the localization region is a 4-volume whose extension 
simply shrinks to zero in physical units, in all directions, in 
the continuum limit. To test which of these two possibilities 
is correct, we define the rms extension Ax of the localization 
region as follows: On a given lattice, let x^ be the site at which 
the density p(x) is largest, and let 



H (A — / ^ol 
W_ \L-|j«f-jtg| 



if < L/2 
otherwise 



Then 



Ax z 



-4- 



-4- 



-4)p( x ) 



(3.3) 



(3.4) 



If the localization region were surface-like, with finite surface 
area in physical units, then unless the surface is space-filling 




FIG. 16: ba 2 in the range (3 £ [2.1,2.4]. If ba 2 is a constant at large 
P, then the localization volume shrinks to zero, in physical units, in 
the continuum limit. 



we expect that Ax is finite in physical units. On the other hand, 
if the eigenmode has its support on a 4-volume of extension Ax 
in all directions (i.e. "point-like"), then we expect that that the 
volume of a 3-sphere of radius Ax (volume = j% 2 Ax 4 ) would 
be comparable to the localization volume b « V /IPR. In fact, 
computing the ratio 



V/IPR 

\k 2 {Ax) 4 



(3.5) 



we find that r| w 0.5 ±0.15, over the whole range of P € 
[2.1,2.4]. We conclude that the localization region is in fact 
point-like rather than surface-like, with a volume that shrinks 
to zero in physical units in the continuum limit. 

The point-like nature of the lowest eigenmodes is confirmed 
by viewing the eigenmode density p (x) for typical thermalized 
lattices. We have made computer animations of p(x), on a 
series of two and three volume slices of the four-dimensional 
lattice, available on the web at the URLs cited in ref. [13] (2D 
slices) and [14] (3D slices). 
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FIG. 17: Plot of lPR/(a z L 4 ) at some individual volumes at higher (3 
values. 



A. The Higgs Phase 

The next question is whether the sharp localization of the 
lowest-lying eigenmode in adjoint representation is due to 
large-scale confining disorder, or whether it is merely due to 
local gaussian fluctuations at the lattice scale. Since the ad- 
joint representation is insensitive to the center of the gauge 
group, it is not possible to cancel out the confinement disor- 
der in this representation by the de Forcrand-D'Elia trick. In- 
stead we can add a Higgs field, and investigate what happens 
to localization as we adjust the Higgs coupling to go from the 
confinement-like to the Higgs phase. For the SU(2) gauge 
group, and the Higgs field in the fundamental representation, 
the action can be written as [15] 

S= p £ ±Tr[UUUW}+yY,^Hx)U,(xMx+!i)} 

plaq x,fj 

(3.6) 

where (j)(x) is an SU(2) group-valued field. It is known that 
(i) this theory has no local, gauge-invariant order parameter 
which can signal a phase transition; (ii) there is no confine- 
ment asymptotically for any y > 0; and (iii) any two points in 
the p-y phase diagram can be joined by a line along which 
the free energy is entirely analytic [16]. There is nonetheless a 
real distinction between the Higgs and the "confinement-like" 
phases. If we fix to Coulomb gauge, then there remains on any 
time-slice a global remnant of the gauge symmetry. This rem- 
nant symmetry is unbroken in the "confinement-like" phase, 
and broken in the Higgs phase [17]. 5 Physically, this implies 
that an isolated charge, together with its Coulomb field, is a 
finite energy state in the Higgs phase, and infinite energy in 
the confinement-like phase. 

At p = 2.1, the transition from the confinement-like to the 
Higgs phase occurs at around y = 0.9. In Fig. 18 we dis- 



play the IPRs at y = 0.7 (confinement-like phase) and y = 1 .2 
(Higgs phase). The drastic reduction in IPR in the Higgs phase 
as compared to the confinement phase is evident. This does 
not mean that in the Higgs phase the magnitude of the lowest 
eigenmode is homogeneous throughout the lattice volume, but 
we do find that the lowest eigenmode becomes substantially 
more extended as the lattice volume increases. The qualita- 
tive difference in lowest eigenmodes of the two phases is clear 
from a plot of the remaining norms at p = 2.1, shown for the 
confinement-like phase (y = 0.7) in Fig. 19, and the Higgs 
phase (y= 1.2) in Fig. 20. 



Gauge-Higgs: p=2.1, Lowest Eigenstate 
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FIG. 18: IPRs in SU(2) gauge-Higgs theory at (3 = 2.1 and y = 0.7 
(confinement-like phase) and y = 1.2 (Higgs phase). 
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FIG. 19: Remaining norms in the confinement-like phase on L 4 lat- 
tices, L = 8 - 18. 

In order to express these differences a little more quantita- 
tively, we note that in an Anderson localized state p (x) decays 
exponentially away from the localization region. This means 
that the sum 



5 The symmetry-breaking order parameter is non-local when expressed as 
a gauge-invariant quantity, so there is no contradiction of the Fradkin- 
Shenker-Osterwalder-Seiler [16] result. 



*p = Lp 1/p « 



(3.7) 



is finite in the infinite volume limit, for any p > 0. In contrast, 
if p(x) has only a power-law falloff (or no falloff at all, for 



8 



Gauge-Higgs, y=1 .2, Lowest Eigenstate 
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FIG. 20: Same as Fig. 19, for the Higgs phase. 



homogeneous states), then there will be some p > 1 such that 
s p — > °o in the infinite volume limit. In Fig. 21 we display 
our results for S2, computed at y = 0.7 and y = 1.2. In the 
confinement-like phase, S2 is constant at approximately S2 — 
71, while in the Higgs phase, S2 appears to rise linearly, on 
L 4 lattices, with extension L. This confirms that the lowest 
eigenmode of the adjoint Laplacian in the Higgs phase is not 
an exponentially localized Anderson-type state. 



Gauge-Higgs Theory, (3=2.1 
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FIG. 21: Test of Anderson (exponential) localization: S2 vs. volume 
in the Higgs and confinement-like phase. 



IV. THE j = 3/2 REPRESENTATION 

We have also calculated the low-lying eigenmode spectrum 
of the covariant Laplacian in the j = 3/2 representation of the 
SU(2) gauge group, at various volumes and p values. We find 
that the low-lying modes are even more strongly localized in 
the j — 3/2 representation than in the adjoint representation, 
and the data falls on a single line when plotting I PR x a 4 vs. 
lattice volume (La) 4 in physical units. This data suggests that 
b itself is approximately p-independent, and if this feature re- 
mains true at higher p, it implies that the localization region 



is always on the scale of the lattice spacing, regardless of the 
length of the lattice spacing in physical units. 

The j = 3/2 representation of SU(2) has the same n-ality 
as the fundamental representation. However, unlike the fun- 
damental representation, vortex removal has essentially no ef- 
fect on localization of the low-lying j = 3/2 eigenmodes, as 
seen by comparing the no-vortex result in Fig. 23 with the 
previous Fig. 22. Given the high degree of localization in 
the j = 3/2 representation, insensitivity to vortex removal is 
not surprising. It is rather obvious that eigenmodes of center- 
projected configurations, in the j = 3/2 representation are in 
complete correspondence to the eigenmodes of the j = 1/2 
representation, and would therefore be localized on the scale 
of fermis. This sort of physical-scale localization, which is re- 
moved in the fundamental representation by vortex removal, 
cannot easily account for the lattice-scale localization found 
in the j = 3/2 representation. 



Cov. Laplacian in j=3/2 rep'n, 1 eigenvector 
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FIG. 22: Log-log plot of IPR x a 4 vs. lattice volume in physical units, 
for the j = 3/2 representation. 
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FIG. 23: Same as Fig. 22, for vortex-removed configurations. 
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V. LOCALIZED MODES VS. MASSLESS PARTICLES 

The propagator of a color-charged scalar particle of bare 
mass OTq, propagating in a background lattice gauge configu- 
ration Ufj(x), can be expressed as 



S ab (x,y;U) 



1 



-A +OTq 



ab 



' K+ml 



(5.1) 



where the (j)JJ(x) are the eigenmodes of of the covariant lat- 
tice Laplacian operator. The operator S ah is gauge covariant, 
and therefore its expectation value vanishes in the absence 
of gauge-fixing, regardless of confinement properties. If a 
physical (e.g. Coulomb) gauge is chosen, then we expect that 
(S ab (x,y;U)) = in a confining theory if xq ^ yo, because a 
color-charged state should not be able to propagate. This re- 
quirement is enforced in Coulomb gauge, as explained in ref. 
[17], by the existence of an unbroken remnant global gauge 
symmetry, which ensures that S ab will still average to zero for 
finite time separations. 

However, suppose we wish to restrict our discussion to 
gauge-invariant operators. Then the next simplest two-point 
function to consider is 



K(x-y) = (S ab (x,y,U)S ba (y,x,U)) 



(5.2) 



which can be regarded as the quenched approximation to the 
VEV of two identical scalar fields cp 1,2 



9 lfl (*)9 t2a (*)9 ie (y)9 T2e (y) 



16/ 



A2b( 



(5.3) 



where the two "flavors" are introduced just to avoid two 
scalars annihilating into an intermediate state of gauge bosons. 

Let us now ask whether the bare mass m 2 , can be adjusted 
so as to make K(x — y) long range. If that is possible, then the 
two scalars must couple to an intermediate state containing 
only massless particles. These massless intermediate particles 
might be 

1. The scalars themselves. In a non-confining theory such 
as QED4, it is generally possible to adjust ml so that the 
physical mass m p h ys of the scalar particles takes on any 
value desired; in particular m pnys = 0. 

2. A two-scalar bound state B, with rag = 0. 

3. "Gluelump" states G. Scalars in the adjoint representa- 
tion can bind to gluons to form colorless bound states. 
In that case, it might be possible to adjust m 2 , so that 
mo = 0. 



We will argue that localization of the low-lying states rules 
out any long-range, power-law falloff (such as l/\x — y\ 4 or 
l/|x — y| 2 ) in K(x — y), at least if tachyonic modes are ex- 
cluded, and hence the mass of the scalar particles cannot be 



adjusted to zero as in non-confining theories. The reasoning 
is closely related to that of McKane and Stone [18] (see also 
Golterman and Shamir [3]), in connection with Anderson lo- 
calization and the absence of massless Goldstone particles. 

The argument goes as follows: We can express K(x — y) as 
a double mode sum 



K(x-y) 



T K(*M b (y) y fe(y)tf (*) 



Km + ml 



(5.4) 



Suppose that we could adjust ml such that m phys or m G were 
massless. Then the two-point function K(x — y) is long-range, 
only falling off like 1 j\x — y\ 4 , so that 



1 



V 



x y 



1 



■L 



ly 

V^(K+m 2 ) 2 
P(X) 



dX 



(X + ml) 



2\2 



(5.5) 



is logarithmically divergent at volume V — > °° (where p(X) 
is the normalized density of eigenvalues). It is assumed that 
there are no tachyonic modes; i.e. ml > —X m i n . Then from 
(5.5) we see that (i) an infrared divergence is possible only if 
Wq = —Xmin, and (ii) if there is a divergence, then it is due en- 
tirely to the low-lying eigenmodes of the covariant Laplacian. 
It is therefore sufficient to restrict our attention to the near- 
zero modes of — A + ml. However, if the near-zero modes are 
localized, then it is impossible for those modes to be respon- 
sible for a long range in K(x — y). Remove those modes, and 
the sum over the remaining modes in eq. (5.5) is finite, so a 
long range can't be coming from the higher modes either. The 
conclusion is that K(x — y) is short-range. Given localization 
of the low-lying eigenmodes, it is just not possible to adjust 
ml such that intermediate states of mass mg, mg, or m p h ys 
become massless, at least not without introducing tachyonic 
modes into the kinetic operator — A + m 2 ,. 

There is still one puzzle. The rhs of eq. (5.5) can certainly 
be infinite if ml = —X m i„, but how is it possible for the lhs 
to diverge if K(x — y) is short range? The answer is that if 
we have Anderson localized states, the infinity in (5.5) is not 
coming from the sum over x,y, because in fact there is a di- 
vergence at any x,y. Notice that the operator inside brackets 
in eq. (5.4) contains a piece (from n — m) 



G(*o0 = £ 



2\2 



(5.6) 



and note that the argument of the sum is positive. Now sup- 
pose first that all states (j)„ are extended, and normalize to 
unity. Then | (j>„ (jc) | <~ 1/VV, which means that the numer- 
ator contributes a factor of 1 /V 2 . This has to be compensated 
by sums over x,y, contributing a factor of V 2 , and then the 
divergence is determined by the integral 



1 



V^(X n +m 2 ) 



2\2 



/l 

u m in 



dX 



p(X) 



(X + m 2 ) 



2\2 



(5.7) 



If we take V — > °° without summing over positions, then 
G(x 7 y) — > 0. So in this case, the divergence (if it exists) is 
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infrared in origin. By contrast, suppose the near-zero modes 
of (— A + OTq) are localized. Presumably the number of these 
modes grows linearly with volume. For each x,y, there is one 
particular state (j)„ which gives the largest contribution to the 
mode sum, and from this contribution 



G(x,y) ■ 



1 exp[ 



\y-z n \)/L] 



b 2 



(X n + m ) 



(5.8) 



where |0„ (x) | is maximized atx = z n , and b is the localization 
volume. Certainly G(x,y) may diverge if X„ + m 2 ) — > 0, but 
this does not imply any divergence in the range of G(x,y) in 
the same limit. 



A. Eigenmodes of the Faddeev-Popov Operator 



In the scenario for confinement advocated by Gribov and 
Zwanziger [19], (K ab (x,y,A)} rises linearly with separation 
at large \x — y\, and this has been confirmed by numerical sim- 
ulation [20]. But this behavior can never come about if M _1 
is short range. It follows, from the previous arguments, that 
if M _1 is long-range then the low-lying eigenmodes of the F- 
P operator must be extended, rather than localized. The F-P 
operator actually has three trivial (constant) zero modes with 
IPR = 1. Fig. 24 shows the IPR vs. volume of the lowest 
non-trivial eigenmode at p = 2.1. The IPR values, even on 
sizable lattices, are rather low, and the growth in IPR with 
volume is very gradual. The low IPRs imply that the lowest 
non-trivial mode is an extended state (at least on the lattice 
sizes shown), despite the similarity of the F-P operator to the 
covariant Laplacian, but in accordance with the reasoning pre- 
sented above. 



An interesting check of the preceding logic is provided by 
the Coulomb gauge Faddeev-Popov (F-P) operator. 6 In the 
continuum, this operator looks superficially much like the co- 
variant Laplacian in D=3 dimensions, i.e. 

M = -V-D(A) (5.9) 

where Dk(A) is the covariant derivative, and A is fixed to 
Coulomb gauge. On the lattice, with SU(2) link variables ex- 
pressed as 

U IJ (x)=b M (x) + ia fl -a (5.10) 
the F-P operator becomes 

K = ^L{M**w +**(*-*)] 
k=i 

- § x . y -M x )-\ x -'My)} 

- ^ e t{VKW-VKW} (5.ii) 

k=\ 

Because of its strong similarity to the covariant Laplacian, 
one might expect that the low-lying eigenmodes of the F-P 
operator would also be localized. However, there is a strong 
argument against this. In Coulomb gauge, the Coulomb en- 
ergy due to a static color charge distribution p a (x) is given 
by 

H cou i = \ jd 3 xd 3 yJ-i~p a (x)jK ab (x,y;A)p b (y)rl (5.12) 
where 

K ab (x,y;A) = [M~ l (-V 2 )M~ l ] J 

j = det[-V-£>(A)] (5.13) 



6 We thank Dan Zwanziger for this suggestion. 



F-P operator, 4th eigenvector 
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FIG. 24: IPR of the lowest non-trivial eigenmode of the Faddeev- 
Popov operator. 



VI. CONCLUSIONS 

The low-lying eigenmodes of kinetic operators on the lat- 
tice, such as the Dirac operator and the covariant Laplacian, 
are probes which may be useful in exploring the properties of 
the QCD vacuum. In this article we have begun a study of 
the low-lying eigenmodes of the Yang-Mills covariant Lapla- 
cian, and have uncovered some intriguing properties of those 
eigenmodes which we cannot as yet explain. 

The most puzzling feature concerns gauge-group represen- 
tation dependence. We have found that eigenmodes of the 
Yang-Mills covariant Laplacian operator are localized at the 
lower and upper ends of the spectrum, and that the degree of 
localization of these eigenmodes depends qualitatively on the 
group representation of the Laplacian operator. With b denot- 
ing the localization volume, in lattice units, of the lowest-lying 
eigenmode, the data indicates that: 

1 . In the fundamental representation, the localization vol- 
ume fe(p)a 4 (p) in physical units is p-independent, with 
an average extension on the order of 2 fm. The local- 
ization property disappears when center vortices are re- 
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moved from the lattice configuration, while the local- 
ization property of center-projected lattice configura- 
tions is similar to that of unmodified configurations. 

2. In the adjoint representation, the quantity Z?(p)a 2 (p) is 
p-independent. This means that at weak couplings the 
localization volume is large in lattice units, but goes to 
zero, in the continuum limit, in physical units. The ad- 
joint representation is insensitive to vortex removal, but 
we find that exponential localization, in a gauge-Higgs 
theory, goes away in the Higgs phase. 

3. In the j — 3/2 representation, b itself is insensitive to 
coupling, which indicates that the localization volume 
is on the scale of the lattice spacing. Vortex removal 
has no effect on this high degree of localization. 

We have also found that localization is absent in the eigen- 
modes of the Faddeev-Popov operator in Coulomb gauge, de- 
spite the close resemblance of the F-P operator to the co- 
variant Laplacian. This absence of localized, low-lying F-P 
eigenmodes is implied by the existence of a confining color- 
Coulomb potential, and is therefore essential in Coulomb 
gauge confinement scenarios. 

Our results cast some doubt on the utility of the so-called 
Laplacian gauges [21], particularly on large lattices. Lapla- 
cian versions of various gauges [22] have been introduced in 
the literature; these have the advantage of avoiding the Gribov 
copy problem. However, each of these gauges makes use of 
either the lowest, or the lowest two or three Laplacian eigen- 
modes, on the grounds that these are uniquely the smoothest 
of the eigenmodes. We have found, however, that the lowest 



eigenmodes are actually the least smooth of all the Laplacian 
eigenmodes, and on a very large lattice are negligibly small on 
most of the lattice volume. Presumably, such highly localized 
modes are not well motivated for the purpose of gauge-fixing. 

The relationship of localization to confinement deserves 
further investigation. In the fundamental representation, 
where localization volume is finite in physical units, there is 
evidence of a strong connection between confinement, center 
vortices and the localization effect. The relationship between 
localization and confinement disorder is not so clear for the 
higher group representations, where the localization lengths 
tend to zero in the continuum limit, and are unaffected by 
vortex removal. It is possible that localized Laplacian eigen- 
modes in different group representations are probing differ- 
ent features, and perhaps different length scales, of the Yang- 
Mills vacuum state. 
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